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We show that observational limits on the possible time variation of constants of Nature are sig- 
nificantly affected by allowing for both space and time variation. Bekenstein's generalisation of 
Maxwell's equations to allow for cosmological variation of a is investigated in a universe containing 
spherically symmetric inhomogeneities. The time variation of a is determined by the local matter 
density and hence limits obtained in high-density geophysical environments are far more constrain- 
ing than those obtained at high redshift. This new feature is expected to be a property of a wide 
class of theories for the variation of constants. 
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There has been renewed interest in modelling and 
constraining possible time variations in the traditional 
constants of Nature. New observational limits have 
been stimulated by high-quality astronomical data [Q, 
whilst theories of high-energy physics featuring addi- 
tional dimensions of space and new dilaton fields have 
provided motivations for studying variations in the gravi- 
tational, strong, and electroweak coupling 'constants' 
1^. There exists a well-defined theoretical framework 
for studying variations in the gravitational coupling by 
means of scalar-tensor gravity theories and these have 
been studied in the case of purely time-dependent varia- 
tions arising in spatially homogeneous cosmological mod- 
els and in the static case. Different theoretical descrip- 
tions admitting variations in the fine structure constant, 
a, have been given by Bekcnstcin and by Barrow and 
Magueijo Although it is usual to refer to all these 
theories as descriptions of a variation in the value of 
a dimensional constant (like G, c, or e), only variations 
in dimensionless constants have an invariant operational 
meaning under changes of units. Thus, it is possible to 
find transformations which map the representation of a 
theory of varying a into either a form with varying c, 
1^, or varying e, j7|.The published observational limits 
on the variations of constants are derived from local geo- 
physical and laboratory considerations , or from high- 
redshift astronomical observations of atomic spectra [Q, 
1^ . Generally, the direct laboratory bounds are sig- 
nificantly weaker than those derived from geophysics and 
astrophysics because of the billions of years of look-back 
time over which the latter two fields can gather data. 
These observations are traditionally used to place limits 
on how much the value of a particular constant of Nature 
can have varied over an interval of time. It is implicitly 
assumed that the rate of time variation is the same ev- 
erywhere. 

In this paper we show why it is important to include 
the possibility of spatial variations in the study of time- 



variations of physical constants. Their inclusion is es- 
sential if local solar system or terrestrial limits on the 
possible variation of constants are to be correctly corre- 
lated to cosmological observations. Spatial variations in 
the values of constants develop by analogy with spatial 
variations in the density, driven by the process of grav- 
itational instability as the universe expands. We show 
how the two sets of variations are related in a partic- 
ular theory which is most expeditious for this purpose. 
This is the lagrangian theory for the variation of the fine 
structure constant developed by Bekenstein |^ . 

Bekenstein has formulated a theory with varying fine 
structure constant, a = e^ /he, defined by the require- 
ments that at constant a it reduces locally to Maxwell's 
equations with minimal vector coupling; all variations of 
a are dynamical and derivable from an action principle; 
the action is causal, time symmetric, and gauge invari- 
ant; no length scale smaller then the Planck length enters 
the theory; and the gravitational field equations are un- 
changed. These requirements lead to an additional con- 
tribution to the action 



where the fundamental electron charge varies through a 
dimensionless universal field £{x'^) defined by e = eQe{x'^) 
with eo constant and £ is a constant length scale defin- 
ing the theory. The variation of this action leads to a 
propagation equation for e : 
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where □ is the covariant D'Alembertian. In the case of a 
spatially homogeneous Friedmann background universe 
of any curvature, with expansion scale factor a(i), this 
simplifies to 
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where dot denotes difFerentiation with respect to the co- 
moving proper cosmological time, t, p is the matter den- 
sity of the universe, and ^ ~ 10~^ is a dimensionless 
measure of the fraction of mass in Coulomb energy for 
an average nucleon. Notice that in the dust-dominated 
era of the universe the right-hand side of (|l|) is constant 
since p cx a~^. From now on we set h = c = 1. The 
variation of the fine structure is Hnked to that of e by 
a/a 2e/e. 

We can use eq. (|l|) to investigate the effects of space 
and time variations in s. Since the background cosmo- 
logical evolution for a{t) is unchanged by introducing 
variations in e we can exploit the Birkhoff property of 
the gravitational sector of the theory to create an exact 
description of a universe possessing spherically symmet- 
ric inhomogeneities in density, curvature and, electron 
charge by matching together solutions of different spa- 
tial curvature and density. We can then calculate the 
evolution of e in both regions and correlate them with 
inhomogeneities in the density. The Friedmann equation 
describing the evolution of a dust inhomogeneity with 
density parameter fip and Hubble parameter Hq is 
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Using 



we 



can integrate (|^), to obtain 
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where we have defined the constant 



where tp is the Planck length. We note that in the flat 
background universe (f2 = 1) we have i3 = 0, so (see also 
0) 



e 3a3/2 



and 



Ine : 
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where a(t) = {t/to)"^/^ and a boundary condition e(to) 
1 is imposed, following B. Hence, we have 



exp[— ^{1- (y)}]. 



In a closed universe, when > 1 we integrate (|^) to 
obtain a generalisation of (^, 
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where C, — — and ci is a constant. 

These results form the basis of a quantitative model 
for the variation of e in time and space. We model an 
inhomogeneous universe containing spherically symmet- 
ric curvature inhomogeneities by matching a flat (k = 0) 
Friedmann universe with scale factor a{t) = R{t) contain- 
ing matter with zero pressure and density p to a curved 
(fc 7^ 0) Friedmann universe with scale factor a{t) = S{t) 
and density p + S The standard metric matching condi- 
tions, R = S and dR/dt = dS/dt, are applied initially 
at t = ti. The hydrostatic support is maintained in the 
absence of pressure and allows the same time coordinate 
to be used in both regions. From the conservation equa- 
tion we have psS^ = (1 -I- 5)pR^ where ps is the density 
in the lump and p is the density in the flat background 
universe and we define S = (ps — p)/p and 5i = S{ti). 
The matching of the fiat and curved regions gives, for 
the overdense region 

= C [{1 + 6i)S-'^ ~ SiS^^] 

Integrating and using the Friedmann equation for the 
background evolution of R{t) cx t^^"^, this gives the rela- 
tion 

^i?3/2(i + ^^)l/2 - S+{S - 1^)1/2 
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where = Si{l+S^^) is the expansion maximum in the 
positively curved lump when JIq > 1, and C is a constant 
that can be expressed in terms of 5*1 and di alone. When 
the density contrast is small {S < 1) we have S cx t^^^. 
Using (|l|), we have the following relationship linking the 
rate of change of the electron charge in the background 
universe, e{t), to its variation in a density inhomogeneity 
with contrast, S, described by Ss{t), 
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Eqn. (H) can also be applied to the case of underdensities 
(Voids') modeled by open Friedmann regions with CI < 1 
when the contrast parameter lies in the range —1<S< 
0. It shows how variations in local cosmological density 
influence are coupled to the time variation of e. The 
effects are likely to be large when comparing intra and 
extragalactic effects because a typical virialised galaxy 
has an overdensity of <5 ~ 10^. The spatial distribution 
of density inhomogeneities can be described by means of 
the 2-point correlation function ^(r) ~ {lMpc/r)^^^ [ p^ . 
Statistically homogeneous and isotropic spatial variations 
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in the statistics of eje should therefore follow the same 
trend with 

< — (x + r) — (x) >cx r~^-^ 

The propagation equation (|^) has the form □(In e) oc p. 
The action of the D'Alembertian on Ine ensures that the 
theory is hyperbolic with nice properties fip^. Typ- 
ically, we expect other theories incorporating varying 
'constants' derived from the variation of a scalar field 
to have a similar structure In the case of Brans- 

Dicke theory we find a similar form, with, the Brans- 
Dicke scalar, ~ G"^ governed by |14[ 
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where T is the trace of the energy momentum tensor 
of the matter sources. Again we see how there is di- 
rect coupling between the local density of matter (on 
the right-hand side of eq. (|^)) and the variations of 
the scalar field associated with G. The full system of 
Brans-Dicke field equations is more complicated to solve 
than those of the Bekenstein theory examined above but 
the same general structure is obtained. A bound over- 
density with large density contrast over the flat back- 
ground universe will possess a different value of 4> and 
of (}) at any given time compared to the values of these 
quantities in the background. In the weak-field limit the 
Brans-Dicke cosmological equations have the decoupled 
linear form for the evolution of the gravitational potential 
perturbation, the density excess over the background 
universe S, and 4" = S(l)/(f> which gives the relative pertur- 
bation of the Brans-Dicke parameter over the background 
0(t), so (5 = when * = 0, 
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As in the Bekenstein theory, the magnitude of the vari- 
ations in G/G from site to site in the universe will be gov- 
erned by the variations in the density. Eqn. (^) shows 
how the density perturbations create inhoniogeneity in 
and hence in G with SG/G S(t)/(f) - $/(2 -I- us). 

Existing observational evidence has assumed that there 
is a single global rate of change in time of a constant, 
like G or a, which is the same everywhere. Theories 
like Bekenstein's for varying a and Brans-Dicke gravity, 
have direct coupling between the variation of the scalar 
whose changes drive variations in the 'constant' and the 
local matter density. The rate of change of the fine struc- 
ture constant in Bekenstein's theory will be proportional 
to the density of the environment where the variation 



is being observed or constrained. From (g) we see that 
the rate of variation is closely linked to the local matter 
density. Thus a constraint on time variation of the fine 
structure constant obtained by examination of the re- 
quirement for a neutron-capture resonance to have been 
in place 1.8 x 10® years ago in the Oklo natural reactor 
under the Earth's surface in Gabon, where p > I 
gmcm~^ corresponds to a cosmological background limit 
that could be 10^® times stronger. Likewise, when one in- 
fers upper limits on the possibility of time- varying G from 
observations in the solar system of the binary-pulsar sys- 
tem they should not be compared directly with cosmolog- 
ical limits without renormalisation for the effects of den- 
sity enhancement due to the formation of galaxies. Thus 
they should correspond to far stronger cosmological lim- 
its when they have been scaled by the appropriate den- 
sity ratio p{solarsystem) / po, where po is the background 
density of the universe. These considerations should pro- 
voke a reassessment of the relationships between different 
astrophysical and geological limits on possible variations 
of constants. The overall result will be to dramatically 
strengthen restrictions on allowed cosmological variations 
of the traditional fundamental constants. Typically, ter- 
restrial and Galactic limits of time variation could be 
strengthened by a factor of order p{local) / po ~ 10^® or 
p{galactic) / Po 10^ when translated onto cosmologi- 
cal scales, depending on the details of the Galaxy and 
planetary formation process. This has significant con- 
sequences for the investigation of additional spatial di- 
mensions that appear to be required in some versions 
of string theory and M theory ||^ and some versions of 
quintessence We hope that this will stimulate a de- 
tailed analysis of the complex sequence of events that 
transform the almost-homogeneous evolution of varying 
constants on an extragalactic cosmological scale into the 
complicated but greatly amplified pattern of variations 
that would be produced by gravitational collapse to the 
high-density islands of matter in which and on which we 
reside and from which we make our observations. 
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